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In classical physics, chaotic behavior constitutes the ma- 
jority of all possible dynamics of isolated conservative sys- 
tems. Because every physical process ultimately follows 
the rules of quantum mechanics, in principle one should 
frequently observe a corresponding quantum chaotic evo- 
lution, with all the signatures of classical chaos £Q. How- 
ever, the fragile quantum coherence effects quickly vanish 
even under weak interactions of the system with its sur- 
roundings, restoring classical behavior. The loss of coher- 
ence occurs even faster for chaotic systems 0, and is hard 
to avoid, even in well-controlled laboratory conditions. As 
a consequence, few experiments have been performed with 
quantum systems that have a classical chaotic counterpart 
I3HI2I. Here we report the implementation of the quan- 
tum Kicked Harmonic Oscillator, a paradigmatic system 
for the study of quantum chaos. Using an all optical setup 
that is essentially decoherence-free, we are able to observe 
the non-linear dynamics through the direct measurement 
of the time-evolved Wigner function. Completely control- 
lable experimental parameters allow us to adjust the sys- 
tem continuously from regular to completely chaotic dy- 
namics, as well as the effective Planck constant, the key 
parameter associated to the quantum-classical transition. 
We show the utility of our approach by experimentally in- 
vestigating the decoherence induced by our chaotic system 
on a quantum bit. Our scheme can be employed in a rich 
variety of experimental investigations of non-linear quan- 
tum systems. 

Chaotic classical systems have the characteristic trait of be- 
ing extremely sensitivity to initial conditions. This behav- 
ior, together with the experimental imprecision of the initial 
conditions, cause these deterministic systems to be inherently 
unpredictable. The question as to how classical chaotic dy- 
namics manifests itself in quantum mechanics is an impor- 
tant branch in the field of quantum chaos. In addition to fun- 
damental questions concerning the correspondence principle 
and the semi-classical limit of quantum mechanics, a number 
of intriguing quantum-dynamical features have been uncov- 
ered. Prominent examples are dynamical localization |[T3lL 
the quantum suppression of classical diffusion, and the en- 
hancement of the tunneling rate in the presence of chaos in 
the corresponding classical dynamics (4). 

The simplest systems that exhibit all the intriguing man- 
ifestations of classical chaos in their quantum dynamics are 
periodic time-dependent Hamiltonian (Floquet) systems, and 
are thus the most vastly studied systems in the discipline of 
quantum chaos fT4l . The quantum evolution up to discrete 



time t = nT is described by the quantum map 

|V(n)> = U n |V(0)>, (1) 

where n is an integer and the Floquet operator U describes 
the unitary quantum evolution in the time period T. They 
have been extensively investigated from a theoretical view- 
point, quantum maps corresponding to classical chaotic be- 
havior have been studied experimentally in different systems 
BUTTIl. This is due primarily to the difficulty in manipulating 
quantum systems and preserving them from decoherence, the 
interaction between a quantum system and the environment 
0. 

There have been a few theoretical proposals to realize quan- 
tum chaotic maps using optics fBlim . To our knowledge, 
these experiments have never been carried out, and more elab- 
orate evolutions may not be possible or require difficult or 
complex optical devices. Here we present an optical imple- 
mentation of the quantum kicked harmonic oscillator (KHO), 
an archetype in the study of quantum chaos. The dynamics are 
observed by direct measurement of the Wigner function, and 
all relevant parameters are easily controlled. We then illus- 
trate the utility of this experimental system studying the loss 
of coherence in different chaotic regimes. 

The quantum KHO serves as a model for charges moving in 
orthogonal magnetic and time-dependent electric fields (T8), 
for the electronic transport in semiconductor lattices lfT9ll20lL 
and for trapped ions addressed by a periodic standing laser 
field I2TI . Its evolution is described by the map ([]]), where the 
iteration operator is 

Ukho = R a Vx. (2) 

The operator R a describes the usual evolution of a quantum 
harmonic oscillator parameterized by a = ujT, where uj is the 
oscillator frequency, and T the interval between the periodic 
perturbations. describes a periodic perturbation that cor- 
responds to a potential if cos (Q + (/>). The quantum KHO 
presents a rich variety of dynamical regimes |22, 23 ], depend- 
ing on the perturbation amplitude, K. For a = 2ir/m, with 
m G {1,2, 3, 4, 6}, the quantum system can present diffusion 
in energy; while for other values of a, an interesting transi- 
tion from dynamical localization to derealization is observed 
in numerical simulations I24ll25l . 

The Poincare section of the classical KHO for a = 2ir/m 
is characterized by the appearance of a "stochastic web" as- 
sociated to the chaotic behavior, with periodic regions corre- 
sponding to essentially regular dynamics in between fT8lL as 
illustrated in Fig. [2] The size of the web and the perturbation 
of the regular dynamics inside the periodic regions is governed 
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by the intensity of the perturbation K. For K < 1 the size of 
the web is considerably small, as is the perturbation of the reg- 
ular dynamics inside the periodic regions. When K = 2 the 
KHO can be considered a weakly chaotic system. 

We implement the quantum KHO in the spatial degrees of 
freedom of monochromatic paraxial light, as an optical ver- 
sion of the operator ([2]). This is based on the isomorphism 
between the paraxial wave equation and the Schrodinger equa- 
tion (26||22]|. The light beam is sent n times through a com- 
bination of optical elements designed to implement the op- 
erator U, using the optical system illustrated in Fig[TJa). In 
this scenario, the transverse position in the near and far-field 
correspond to the position and momentum of the photons in 
the beam, and are analogous to the transverse position Q and 
transverse momentum P of a quantum particle. The instanta- 
neous "kick" perturbation, 

K = e h ea , (3) 

is produced using a high-definition Holoeye spatial light mod- 
ulator (SLM). We define Q = vq as the dimensionless version 
of the near-field variable q. v is the spatial period of the cosine 
function, K is the kick strength, and the effective Planck con- 
stant fi e ff will be defined below. The SLM is a relatively new 
device, consisting of a high-definition LCD screen that can 
be used to imprint a phase exp[z/(x, y)] on an optical beam, 
where f(x, y) is a programmable function. As such, the SLM 
allows for the implementation of a number of interesting dy- 
namical maps. 

The harmonic evolution operator R a produces a phase 
space rotation by an angle a = ujT. This is equivalent to a 
Fractional Fourier Transform (FRFT) of order a [28], 

R a = e W2 e ^«(Q 2 +r 2 )/2. (4) 

The dimensionless momentum variable is P = vf a, where 
a is the spatial frequency and f = / sin a is the fractional 
focal length. The optical FRFT system consists of a lens of 
focal length / placed between two sections of free- space of 
distance z a = 2f sin 2 (a/2). The dimensionless operators Q 
an P obey the commutation relation l26ll27lL 

[Q,P] =^eff, (5) 

where H e ^ = v 2 f'/kis the dimensionless effective Planck 
constant, and k is the wavenumber. By changing the order a 
of the FRFT (harmonic evolution between kicks), the ampli- 
tude K of the cosine phase implemented with the SLM (kick 
amplitude), and the spatial frequency v of this phase (effective 
Planck constant), one can manipulate all the relevant parame- 
ters in the dynamics of the KHO. 

The complete experimental setup is illustrated in Fig[TJb). 
A 632. 8nm He-Ne laser is coupled to a single-mode optical 
fiber to obtain a Gaussian beam at the input of the KHO setup. 
This defines a Gaussian light beam as the initial state 1^(0)) 
which is then evolved by the quantum KHO propagator ^ in 
one dimension of the transverse spatial degrees of freedom. 
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FIG. 1 : Optical scheme and experimental setup, a, illustration of 
an optical system to implement the quantum KHO. b, experimental 
setup for the implementation of the quantum kicked harmonic oscil- 
lator and measurement of the spatial Wigner function. 

The state is sent through n iterations of the KHO operator 
Ukho by reflecting n times between the SLM and a mir- 
ror. A polarizing beam splitter (PBS1) is used to polarize the 
beam parallel to the active axis of the SLM display. The har- 
monic evolution between kicks corresponds to propagation in 
the a-order FFT system, which consists of free space propa- 
gation and the cylindrical lens. For practical reasons, we ac- 
tually implement two consecutive FFTs of order a/ 2, before 
each incidence on the SLM. The focal length of the cylindrical 
lens is / = 150mm, and the free space propagation length is 
z = 75mm, so that a = tt/3 between two consecutive kicks 
(see Fig[TJb)). Since the entire KHO evolution is made with 
optical elements that act only in one spatial dimension, the 
perpendicular direction evolves according to free space prop- 
agation. 

To characterize the chaotic dynamics, we do a point-by- 
point direct measurement of the optical Wigner function (28J 

W(Q,P) = aJ°° * (q + (q - e ih ^d^ 

(6) 
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FIG. 2: Experimental and theoretical Wigner distributions for a single iteration of the KHO map and associated classical Poincare 
sections, a-c, directly measured optical Wigner distributions, e-g, corresponding numerical predictions. For (a) and (e), the relevant parameters 
are K — 7.4 and h e ft = 4.72. For (b) and (f), K — 7.4 and h e ft = 0.9. The yellow line outlines the phase space manifold that is the skeleton 
of both the quantum and the classical distributions corresponding to the KHO Hamiltonian, c, g, For K = 2 and /i e ff = 0.9, d, h, the associated 
classical phase space is illustrated by the Poincare sections. In all cases the harmonic evolution parameter is a = tt/3. 



where A is a normalization constant. We use an interferomet- 
ric method l30l to measure the Wigner function, as described 
in the supplementary material. 

Two spherical lenses (LI and L2) of focal length / = 350 
mm are used to map the output state of the KHO system 
\ip(n)) from the transverse plane at position z to the input 
mirror of the Sagnac interferometer. A Dove prism (DPI) 
tilted at a 45° angle is used to swap the horizontal and verti- 
cal coordinates, because for convenience the quantum kicked 
Hamiltonian is implemented in the vertical axis, while the 
Wigner measurement is performed in the horizontal. The sec- 
ond PBS (PBS2) and the first half wave plate (HWP1) are used 
to keep the beam linearly polarized before entering the Sagnac 
interferometer. 

Figs. [2ja-c) show three experimental Wigner distributions 
for a single iteration (n = 1) of the quantum map ^ ap- 
plied to an initial squeezed Gaussian state |^(0)) centered at 
the phase space origin, which is a fixed point of the classical 
KHO map. The theoretical Wigner distributions are shown 
for comparison in Figs. [2|e-g). All three cases correspond 
to the harmonic evolution a = tt/3. Figs. [2ja, e) and (b, f) 
show results for the kick amplitude K = 7.4. One can see 
that the smaller h e ^ = 0.9 used in Fig. [2jb, f) enables finer 
interference patterns to appear, in comparison to Fig. [2ja, 
e) with fi, e ff = 4.72. In the chaotic regime, stretching and 
folding of the Wigner distribution W(Q, P) gives rise to fine 
oscillatory phase-space structure [31 ], known as sub-Planck 
structure, in a scale that is proportional to f^ ff . These in- 



terference patterns are a quantum signature that can only be 
wiped out by decoherence (32j[33]|. However, as the effective 
Planck constant approaches the semiclassical limit (ft e ff —> 0), 
one identifies the concentration of the phase-space distribu- 
tion over the evolved classical curve (yellow line), that is the 
skeleton of both the Wigner function and the corresponding 
classical probability distribution under KHO evolution of an 
initial Gaussian wave packet I34ll35ll . 

The quantum sub-Planck and underlying classical struc- 
ture associated to the evolution of Gaussian wave packets in 
chaotic systems can be observed for two final states corre- 
sponding to weak chaos K = 2 in Figs. |2jc, g) and strong 
chaos K = 7.4 in Figs. |2jb, f). We notice the typical stretch- 
ing and folding of the classical manifold that supports the 
Wigner distribution in these chaotic regimes J34} [35]|. Every 
two points on this classical manifold generate an interference 
fringe halfway between the chord that joins them with a fre- 
quency that is proportional to the length of this chord l34l . By 
comparing Fig|2jb,f ) to Fig|2jc, g), one can observe that the 
skeleton of the Wigner function after one kick for K = 7.4 
folds in phase space in a more pronounced manner due to the 
more chaotic dynamics (see the corresponding Poincare sec- 
tions in Fig. [2jd) and[2jg)). 

Fig. [3] displays experimental Wigner distributions for the 
first three iterations of the KHO map ^ for the same initial 
state |^(0)). In this case we changed the harmonic evolution 
to a = 27r/3. Figs. [5Ja-c) show the evolution for essen- 
tially regular dynamics, with ft e ff = 0.42 and K = 0.75. The 
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FIG. 3: Experimentally obtained Wigner distributions corresponding to the first three iterations of the KHO map. a-c, The case of 
regular dynamics, with K — 0.75, and ft e ff = 0.42. d-f, Weak chaos regime, with K — 2, and h e s = 0.9. In all cases a = 2tt/3. 



Wigner function, initially centered at the origin, spreads over 
a limited area in phase space. Conversely, dynamical spread- 
ing of the Wigner distribution is observed for the weak chaotic 
regime in Figs. [3jd-f) with ft e ff = 0.9 and K = 2. 

As an example of the utility and versatility of our approach, 
we experimentally studied the loss of coherence by a polar- 
ization qubit interacting with a chaotic environment. The re- 
sults are shown in Fig. [4] (see supplementary material for 
methods). Taking advantage of the fact that the SLM only 
imprints a phase on the horizontal polarization component, a 
dephasing-type coupling was implemented between the po- 
larization and the transverse spatial degrees of freedom of the 
laser beam. Decoherence produced by the interaction with 
only a few environmental degrees of freedom is the main ob- 
stacle to quantum computing when the system of interest is 
well isolated 1 36 -38]. In this case, entropy production and ir- 
reversible loss of coherence are directly related to chaotic dy- 
namics of the environmental degrees of freedom |[38l - l4T1l . 

In this case, the quantum KHO, implemented in the spa- 
tial degrees of freedom acts as an environment for the qubit 
encoded in the polarization state. By preparing the incoming 
beam in the linear diagonal polarization state, the polarization 
and the transverse spatial variables become entangled, lead- 
ing to loss of purity of the polarization state. We studied this 
process as a function of K and ft e ff • 

In Figs. |4ja) and (b), the purity of the polarization state is 
plotted as a function of the number of iterations of the KHO 
map. The different curves correspond to different values of 
ft e ff. Fig. |4ja) shows the case in which the KHO has reg- 
ular dynamics (K = 0.5), and Fig. [4] (b) is for essentially 
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FIG. 4: Purity loss of the diagonal polarization state for regular 
and chaotic dynamics of the environment. a,b, Purity as a function 
of the number of KHO iterations, n, for regular dynamics (K = 
0.5), and chaotic dynamics (K = 2), respectively. In both figures 
a = 27r/3, and ft e ff = 0.05 (green triangles), 0.1 (red squares), 0.5 
(purple crosses), 1.0 (orange circles), and 1.5 (blue polygons), c, 
The final purity as a function of ft e ff for n — 3 KHO iterations, with 
K — 2 (green triangles) and K — 0.5 (red squares). 



chaotic dynamics (K = 2). In the latter, for all values of 
h e fi , rapid loss of purity occurs while the state becomes nearly 
maximally entangled (purity equal to 1/2). Fig. |4jc) shows 
the final purity for n = 3 kicks as a function of h e ^ for if = 2 
(green triangles) and K — 0.5 (red squares). The results show 
that the entanglement between the qubit and the environment 
depends strongly on the effective Planck constant for regular 
dynamics (K = 0.5), and this dependence is much weaker for 
the chaotic case (K = 2). 

Our approach allows for the study of the dynamics of a non- 
relativistic quantum system using an intense classical laser 
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beam due to the analogy between quantum mechanics and 
classical wave mechanics. A possible next step is to study the 
chaotic evolution of entangled photons. The optical realiza- 
tion of non-linear quantum dynamics should prove invaluable 
in the experimental investigation of quantum chaos, decoher- 
ence, and the quantum-classical boundary. 

Supplementary Material 
Method to Measure the Optical Wigner Function 

The method used to directly measure the optical Wigner 
function is an interferometric scheme proposed in reference 
(30). The interferometer is illustrated in Figure 1 of the 
manuscript. The displacement and tilting of a steering mirror 
(Ml) at the entrance of a three-mirror Sagnac interferometer 
displaces the optical field by Q and changes its direction of 
propagation by P (which in the paraxial approximation cor- 
responds to the addition of a phase). This produces the field 
&v(ik eS P£/2)V(Q + £/2,z). 

A polarizing beam-splitter divides the field into two spa- 
tially identical components and a Dove prism (DP2) placed 
inside the interferometer realizes opposite 90° spatial rota- 
tions in the two counter-propagating transverse spatial modes, 
resulting in a total relative rotation by 180°. The modes are 
recombined, and projected onto the diagonal polarization di- 
rection before detection by an area-integrating "bucket" de- 
tector. The measured intensity is composed of three terms 
/ = Ii + I2 + /int. where the sum I\ + I 2 is constant and 
equals half the total input intensity. The term ii nt originates 
from the interference between the counter-propagating beams. 
Due to the relative rotation implemented by DP2, ii nt is pro- 
portional to: 

W(Q, P ) = A J * (Q + I ) ^* (Q ~ \ ) e ih °« p *dt, 

(7) 

The right-hand side of Eq.(|7]) is an integral over spa- 
tial variable £ of the overlap between the displaced field, 
exp(zS e ffP£/2)^ r (Q + £/2,z), and its complex conjugate 
with the transformation £ — » — £. This transformation cor- 
responds to the 180° relative rotation implemented by DP2. 
It is important to note that astigmatic effects introduced by the 
Dove prism l42ll inside the Sagnac interferometer does not al- 
ter the measured Wigner function l43l . 

In this way, the amplitude of the optical Wigner function at 
point (Q, P) can be obtained by measuring the intensity at the 
interferometer output, , for different settings of the tilt angle 
and displacement of the steering mirror Ml. In our exper- 
iment these parameters were controlled with high resolution 
motorized stages. At the exit of the interferometer we use a 
quarter wave plate (QWP) that is tilted to correct the polar- 
ization aberrations introduced by DP2 inside the interferom- 
eter l44l . Detailed calculations of the light propagation in the 
interferometer, and how the optical spatial Wigner function 



arises from the area integrated intensity measurements can be 
found in El. 

Analysis of the Experimental Data 

The final state of the optical KHO is obtained at the output 
plane zo, indicated in Figure 1(b) of the manuscript. How- 
ever, the optical wave-function propagates through free space 
and several linear optical elements before reaching the steer- 
ing mirror at the entrance of the Sagnac interferometer, where 
the optical Wigner function is measured. One must therefore 
take into account this evolution before comparing measure- 
ments with the theoretical Wigner functions, obtained from 
numerical calculation of the quantum KHO evolution. This is 
done by calculating the total linear transformation M, result- 
ing from propagation through all of the optical elements and 
free space between the output plane zo and the steering mirror 
shown in Figure 1 (b). 

The expected Wigner function in the measurement plane 
can be written as 

W(x,p x ) = W x KHO \x',p' x ) x W^ uss \y',p' y ), (8) 

where (V y' p' x p' y ) T = M _1 (x p x 0) T are the trans- 
formed coordinates. The function W x KB ~°\x,p x ) represents 
the Wigner function due to KHO evolution of the initial Gaus- 
sian state implemented in the x transverse spatial direction. 
The function Wy Gauss \y , p y ) refers to the Gaussian state de- 
scribing the y transverse spatial direction at plane zq. This 
spatial mode does not undergo KHO the evolution. 

Experimental errors associated to the propagation though 
the optical elements between the output plane zo and the steer- 
ing mirror results in scaling, skewing and rotation of the fi- 
nal measured Wigner Function, W^ exp \x,p x ), in compari- 
son to W( KUO \xip x ). These uncertainties are due princi- 
pally to errors in lens placement, misalignment of the optical 
elements, and diffraction. Nevertheless, a single linear trans- 
formation E corrects these errors, such that W^\x,p x ) = 
W (kiio) (e^^)). Both M and E depend only on the ex- 
perimental setup, and are the same for all measurements of 
Wigner functions, including the case in which the spatial light 
modulator is turned off. In this case, the implemented evo- 
lution is that of a simple harmonic oscillator (SHO). With the 
KHO turned off, we determined the value of E, and used these 
values to correct all of the KHO Wigner functions. This was 
repeated for each harmonic evolution used. 

Using the spatial light modulator to implement a 
dephasing-type decoherence channel 

The spatial light modulator (SLM) imprints a phase on 
the horizontal polarization component of the light beam, but 
leaves the phase of the vertical polarization component un- 
changed. Without the SLM, the optical system is designed to 
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implement SHO evolution via the fractional Fourier transform 
systems composed of the cylindrical lens and free- space prop- 
agation. Therefore, if diagonally polarized light is used in the 
KHO optical setup, one obtains the transformation: 

|*>®|+> ^(Ukho|^)0|^) + U S ho|^)^|^))(9) 

where |+) = (\H) + \ V))/y/2, and \Sfr) designates the trans- 
verse spatial mode of the light beam. The evolution opera- 
tors Ukho and Usho act on the spatial degree of freedom, 
and correspond to the evolution of the KHO and the SHO, 
respectively. The total evolution given in Eq.([9]) could be in- 
terpreted as the quantum evolution of a qubit interacting with 
a quantum KHO via a dephasing-type of coupling of the form 
(KcoaQ)/2\H)(H\. 

Performing a polarization measurement of the output beam 
using a "bucket" detector is equivalent to tracing over the spa- 
tial degree of freedom, and yields, 

p = \ (\H){H\ + \V){V\ + f\H)(V\ + r\V)(H\) , (10) 

where / = (^|Ug HO UKHo|^) is the overlap between the 
SHO and KHO quantum states. The purity of the polariza- 
tion state is given by Trp 2 = (1 + |/| 2 )/2, and decays with / 

EH. 

We use wave plates and a polarizing beam splitter, to per- 
form quantum polarization state tomography using the stan- 
dard recipe [ 45 , 46 ] to obtain the density matrix p, from which 
we calculate the results shown in Figure 4. 
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